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Abstract. Hyperbolic-parabolic systems have spatially homogenous stationary states. When 
the dissipation is weak, one can derive weakly nonlinear-dissipative approximations that gov- 
ern perturbations of these constant states. These approximations are quadratically nonlinear. 
When the original system has an entropy, the approximation is formally dissipative in a natural 
Hilbert space. We show that when the approximation is strictly dissipative it has global weak 
solutions for all initial data in that Hilbert space. We also prove a weak-strong uniqueness 
theorem for it. In addition, we give a Kawashima type criterion for this approximation to be 
strictly dissipative. We apply the theory to the compressible Navier-Stokes system. 



1. Introduction 

We consider hyperbolic-parabolic systems over the 27r-periodic domain T'^ that have the form 

(1.1) dtU + V,-F{U) = V.-[D{U)-VM]. 

These have spatially homogenous stationary solutions. When the dissipation is weak, one can 
derive a weakly nonlinear-dissipative approximation that governs perturbations U about any 
constant solution Uo- These approximations have the form 

(1.2) dtU + AU+QiU,U)=VU, 

where A = FuiJJo) ■ Vx is the linearization of the convection operator about Uo, while quadratic 
operator Q and the linear operator T) are formally given by 

Q(F,F) = hm ^ f e'^Vx- [\Fuu{Uo){e''^Y,e-'^Y)\ dt , 

(1.3) -^-^ 



1 



T 



VY = lim / e'^D{Uo):W^{e~'^Y)dt 



T 



Such weakly nonlinear-dissipative approximations arise when studying incompressible limits of 
the compressible Navier-Stokes system [SI [32], global regularity of fast rotating Navier-Stokes 
and Euler equations [U, [2, [3], asymptotic limits in equations of geophysical fluid dynamics 
[TTj [12] . and fast singular limits of hyperbolic and parabolic PDE's [351 (HI 115] . 

We show that if the original system (11.11) has a thrice different iable convex entropy structure 
then the approximating system (11.21) is formally dissipative in the Hilbert space EI whose inner 
product is given by 
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where HuuiUo) is the Hessian of the strictly convex entropy density H{U) at Uo- This dissipa- 
tion property follows because the entropy structure implies that A is skew-adjoint in H, that 
T> is nonpositive definite in H, and that Q formally satisfies the cyclic identity 

(1-4) = [U I W))^+{V I Q{W, U))^+{W I Q{U, V))^. 

We show that if is also strictly dissipative then the approximating system (11. 2p has a Leray- 
type global weak solution for all initial data in H. We cannot establish the uniqueness of 
these solutions. Indeed, when (11.11) is the Navier-Stokes system of gas dynamics then (11.21) 
includes the incompressible Navier-Stokes system as a subsystem. The uniqueness question 
therefore includes the uniqueness question for the incompressible Navier-Stokes system. We 
can however use the cyclic identity (11.41) to prove a so-called weak-strong uniqueness theorem 
for the approximating system (II. 2p . 

It will be easily seen that T) will be nonnegative definite if and only if the linear operators 
A = FuiJJo) ■ Vx and V = D{Uo) ■ satisfy the Kawashima condition: 

(1.5) no nonconstant eigenfunction of A is in the null space of T>. 

In one spatial dimension it is known that the Kawashima condition implies that V is strictly 
dissipative [25]. We give a stronger Kawashima-type criterion for V to be strictly dissipative 
in higher dimensions. 

Our paper is laid out as follows. Section 2 presents the entropy structure we will impose on 
system (11.11) . Section 3 presents the weakly nonlinear-dissipative approximation (II. 2p and its 
relation to the Kawashima condition (II. 5p . Section 4 develops the properties of the averaged 
operators (II. 3p that we will need later. This includes a proof of the cyclic identity (II. 4p . Section 
5 presents our Kawashima-type criterion for V to be strictly dissipative. Section 6 contains our 
existence and weak-strong uniqueness theorems. Finally, section 7 applies the theory to the 
compressible Navier-Stokes system of gas dynamics. 

2. Hyperbolic-Parabolic Systems with Entropy 
We consider hyperbolic-parabolic systems over T'^ in the divergence form 

(2.1) dtU + Vx- F{U) = V,- [D{U) ■ VxU] , 

where U{x, t) is a density vector over (x, t) eJ'^ X M+ that takes values in C M^. Here 
qpd _ jjd^ (27rZ)'^ is the 27r-periodic torus and is the closure of a convex domain U CM^. We 
assume that the flux tensor F : U ^ W^^^ is twice continuously differentiable such that 

(2.2) dtU + V^- F{U) = is hyperbolic , 

while the diffusion tensor D -.U —>■ ]]j'^x'ixNxN jg continuously differentiable such that 

(2.3) dtU = V^- [D{U) ■ VxU] is parabolic . 

Recall that system (12. 2p is said to be hyperbolic if for every U E U and every ^ G M'^ the 
NxN matrix Fu{U) is diagonalizable within the reals — i.e. it has a complete set of real 
eigenvectors. System (12.20 is said to be strictly hyperbolic if moreover the eigenvalues of 
Fu{U) are distinct. Recall that system (12.30 is said to be parabolic if for every U eU and 
every ^ e M*^ the NxN matrix D(U) : is diagonalizable within the reals and has nonnegative 
eigenvalues. System (12.30 is said to be strictly parabolic if moreover the eigenvalues of D{U) : 
are positive. Many studies of hyperbolic-parabolic systems assume that system (12.20 is strictly 
hyperbolic while system (12.30 is strictly parabolic. We will not do that here. Rather, we 
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will assume that system (12. ip has a strictly convex entropy and satisfies certain nonsingularity 
conditions. 

2.1. Entropy Structure. We say that if : W ^ R is a strictly convex entropy for the system 
(12.11) when H is twice continuously differentiable over U and for every U eU 

(i) HuuiU) is positive definite , 

(2.4) (n) Huu{U) Fu{U) ■ ^ is symmetric for every ^ G M'^ , 

(iii) Hjjjj{U) D{U) is symmetric and nonnegative definite . 

The existence of such a strictly convex entropy implies that system fl2.1l) is hyperbolic-parabolic. 
The compressible Navier-Stokes system is a hyperbolic-parabolic system that is neither strictly 
hyperbolic nor strictly parabolic, yet has a strictly convex entropy. We will study this example 
in Section [71 There are many other systems from physics that fit into this framework [5]. 

Condition (ii) in (12. 4p implies there exists J : U ^ that is twice continuously differentiable 
such that 

(2.5) Hu{U)Fu{U) = MU) . 

It follows that if U {x, t) is a classical solution of (12. ip that takes its values in U then it satisfies 

(2.6) dtH{U) + V,- J{U) = V, ■ [Hu{U) D{U) ■ V,U] - V,Hu{U) ■ D{U) ■ V^U . 
Condition (iii) in (12. 4p implies that for every differentiable f/ : T'^ ^ with U{x) G U 

(2.7) V.HuiU) ■ D{U) ■ V.f/ = V.f/^ ■ HuuiU)D{U) ■ V.f/ > . 

We thereby see that (12. 6p is a local dissipation law for H{U). When (12. 6p is integrated over T*^ 
one obtains the global dissipation law 



(2.8) 



^ j H{U) dx = -J V^Hu{U)-D{U)-V,Udx 

= - J V.U^ ■ Huu{U)D{U) ■ V^Udx < 



2.2. Nonsingularity Condition. The class of stationary classical solutions of (12. ip that take 
values in U is constrained by the entropy structure. It follows from (12. 8p that every such 
solution satisfies 



J V.f/^ • Huu{U)D{U) ■ V,Udx = . 



which by (12. 7p implies that 

V^U"" ■Huu{U)D{U)-V,U = Q. 

Conditions (i) and (iii) then imply that D{U) - ^xU = 0, which when plugged into (12. ip with 
dtU = yields Vx - F{U) = 0. We thereby see that every classical solution of (12. ip that takes 
its values in U satisfies 

(2.9) Vx-F{U) = 0, D{U)-VxU = 0. 

In general this is not enough information to conclude that [/ is a constant. 
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We consider hyperbolic-parabolic systems over T'^ that are nonsine;ular in the sense that for 
every U & U and every continuously differentiable [/ : T'^ — one has that 

(2.10) ^ ' \^ \ ^ Va:t/ = 0. 

D{U) ■ V.f/ = J 

With this additional assumption we can prove the following. 

Lemma 2.1. // the hyperbolic-parabolic system ^2.1\) has a strictly convex entropy (2^_) and 
satisfies the nonsingularity condition ^2.1Ch) then its only stationary, spatially periodic classical 
solutions that take their values in U are constant solutions. 

In the next section we will strengthen the nonsingularity condition f l2.10p . 

3. Weakly Nonlinear-Dissipative Approximations 

We now consider weakly nonlinear- dissipative (WND) approximations to hyperbolic-parabolic 
systems of the form fl2.ip . In order to see how these approximations depend on the choice of 
dependent variables, we express (12. ip in the form 

(3.1) dtU{W) + V,.- F{U{W)) = V,- [D{U{W)) ■ v,.f/(iy)] . 

Here the mapping U : W ^ W is assumed to be a twice continuously differentiable bijection 
and have a nonsingular Jacobian. This means that its inverse mapping U^^ : U ^ W will also 
be a twice continuously differentiable bijection and have a nonsingular Jacobian. 

3.1. Linearization and the Kawashima Condition. In order to motivate our approxima- 
tion, let us first consider the linearization of the hyperbolic-parabolic system (13. ip about some 
constant state Wo G W. This is 

(3.2) dtW + Ao ■ V.W = Bo : V^W , 
where Ao G R'^xNxn^ ^ ^dxdxNxN defined by 

(3.3) Ao = R^'Fu{Uo)Ro , Bo = R-'D{Uo)Ro , 
with UoeU and i?^ e M^'^^ given by 

(3.4) Uo = UiWo) , Ro = dvUiWo) . 

If : W — » M is a strictly convex entropy for system (12.1 p then system (13. 2p is symmetrized 
by the positive definite matrix 

(3.5) Go = RlHuuiUo)Ro. 
Specifically, one sees from (I2.4p and (13. 3p that 

GoAo ■ i = RoHuu{Uo)Fu{Uo) ■ ^Ro is symmetric for every ^ G M"^ , 

GoBo = R^Huu(Uo)D{Uo)Ro is symmetric and nonnegative definite . 

The solutions of (13. 2p thereby satisfy the local dissipation law 

dt{^^GoW) + v. - {\W^GoAoW) = v.- (W^GoBo ■ W.W) - W.W'^ ■ GoBo ■ W^W . 
When this equation is integrated over T'^ one obtains the global dissipation law 

^ / W^GoW dx = - I VrW^ ■ GoBo ■ WrW dx<0. 



dt 
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The initial-value problem for system (13.21) is therefore naturally well-posed in the Hilbert space 
EI = L2(dx; M^) equipped with the inner product 

where Go is given by (13.51) . 

By arguing as we did in the nonlinear settings, every stationary, spatially periodic classical 
solution of system (13.21) must satisfy 

Ao ■ Vjw = , Bo- V^W = . 

We can conclude that X^VT = if system (12. ip satisfies the nonsingularity condition (I2.10p . 
In that case the only stationary, spatially periodic classical solutions of system (13.21) are the 
constant solutions. 

It is natural to ask if every solution of system (13. 2p in HI will decay to a constant solution as 
t — > oo. Kawashima [23l [23] gave an elegant characterization of when this is the case in terms 
of the skew-adjoint operator A and the self-adjoint operator V that are formally given by 

(3.6) ^ = A„ ■ V, , V = Bo:V^. 

Because .4 = ■ \^ is skew-adjoint in H, its spectrum, Sp(^), is purely imaginary. When the 
spatial domain is bounded Sp(^) only contains eigenvalues. Moreover, these eigenvalues and 
their corresponding eigenfunctions are easily computed by Fourier methods when the spatial 
domain is T"^. Kawashima gave the following characterization. 

Theorem 3.1. Every solution of system li^.^) in EI will decay to a constant solution as t ^ oo 
if and only if 

(3.7) no nonconstant eigenfunction of A is in the null space ofV. 

Remark. The Kawashima condition (13.71) is clearly necessary for every solution of system (13.20 
in EI to decay to a constant solution as t oo. Indeed, if V is nonconstant eigenfunction of A 
for the eigenvalue iu such that V lies in the null space of T> then the real part of e~*'^*\^ is a 
real- valued solution of (13. 2p that does not decay to a constant solution as t — >• oo. 
Remark. The Kawashima condition (13. 7p is stronger than our nonsingular condition (I2.10p . It 
has been used to obtain similar results regarding the existence, regularity, and decay as t — >• cxd 
of global solutions to nonlinear systems |3l [ISl [ISl ESI ESI EH EZ] ■ Villani has generahzed it 
in his theory of hypercoercivity [38j where he discusses its relation to the older Hormander 
hypoellipticity condition ^20j|. 

3.2. Weak Dissipation, Two-Time Asymptotics. It is more interesting to consider regimes 
in which the dissipation is weak. Introduce the nondimensional (small) parameter e > so that 

(3.8) dtU{W,) + F{U{W,)) = eV.- [D{U{W,)) -V.UiW,)] , 

Let WoeW and set ^ ^ 

W, = Wo + eW + e^VT^^) + q(^^3^ _ 

Then 

U{W,) = Uo + eRoW + e^Ro{W^^^ + So{W,W)) + 0{e^) , 
where Uo and Ro are given by (13. 4p while So is defined by 

So{W,W) = \Ro^dwwU{Wo){W,W). 
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To first order in e we see that W satisfies 

(3.9) dtW + Ao-V,W = 0. 
The solution of this equation is given by 

W{t) = e-'^W''' , where ^ = A„ ■ . 

Because A is skew-adjoint on H, the solution operator e"*-^ is strongly continuous, one pa- 
rameter group of unitary operators on H. This approximation cannot be valid uniformly in 
time because (1) the solutions of (13. 9p do not decay as t — oo and (2) there are generally 
nonconstant stationary solutions of (13.91) when d > 2. 

In order to overcome these problems, one has to introduce a slow time scale t = et into the 
asymptotics and consider 

edrU{W,) + dtU{W,) + F{U{W,)) = eV.- [D{U{W,)) -V.UiW,)] . 
To first order we see that W still satisfies (13.81) . Hence, 

(3.10) W^(t, r) = e-*^F(r) , 

where the r dependence of Y has yet to be determined. 
To the second order in e we see that 

(3.11) dtW^^'> + AW^^'> = -drW - V^- Qo{W, W) + Bo: V^W - {dt + A)So{W, W) , 
where 

(3.12) Qo{W,W) = \R-'Fuu{Uo){RoW,RoW) . 

The right-hand side of (13. lip is an almost periodic function of t. For W^^^ to be an almost 
periodic, we must require 

lim / e*-^ [drW + V. ■ Qo{W, W)-Bo: V,' W^] dt = . 

Hence, because W is given by (13.100 . we see that Y{t) satisfies 

(3.13) drY +QiY ,Y) =VY , 
where the operators Q and V are formally defined by 

Q{Y,Y)= lim r e'^Q{e-'^Y,e-'^Y)dt, 
(3-14) _ 

VY = lim — / e^^Ve-^^Y dt , 



T 



T^oo 2T „ 

with the operators Q and V being given by 

Q(W, W) = V,.- Qo{W, W), VW = B,: V^W . 
It is easily checked from formulas (13.141) that 
(3. 15) Q{e-'-^W, e-'^W) = e'^^QiW, W) , Pe"*-^ = e~'^~^ 

It then follows from fIXTUD and flXT^ that W{t,et) satisfies 

dtW + A^ + eQ{W, W) = eVW . 
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Setting e = 1, we call this the weakly nonlinear-dissipative approximation of the hyperbolic- 
parabolic system (13. ip . 

3.3. Change of Dependent Variables. Suppose we had expressed the system (13. 6p in terms 
of different dependent variables W: 

(3.16) dtU'iW:) + V,- F{U'{W:)) = e V,.- [D{U'{W:)) ■ V,U'{W^)] , 

where U' : W — > W is a twice continuously differentiable bijection with a nonsingular Jacobian. 
Let be the unique constant state such that U'{W^) = Uo = U(Wo)- If we approximate 
solutions of (I3.16P that are near just as we approximated solutions of (13.80 that are near 
Wo then to leading order we obtain 

dtW' + A',-V,W' = 0, 

where A' = T-^A^T^ with 



(3.17) To = dw'[U-\U'{W'))) 



W'=W' 



In other words, the leading order approximation of the transformed system (13.160 is the leading 
order approximation of the original system transformed by the associated linear change of 
variables W = ToW'. 

Remarkably, the same transformation property holds for the weakly nonlinear-dissipative 
approximation! Specifically, if the weakly nonlinear-dissipative approximation of the primed 
system (I3.16P is 

dtW' + A'W' + eQ'{W', W') = eV'w , 
then it is related to the weakly nonlinear-dissipative approximation of the unprimed system by 

(3.18) A' = T;'ATo, Q\w\W')=T-'Q{ToW\ToW'), v' = T~'VTo, 

where To is again given by (I3.17p . In other words, it is simply the original nonlinear-dissipative 
approximation transformed by the associated linear change of variables W = TqW. This fact 
allows us to derive the weakly nonlinear-dissipative approximation of system (13. 8p in any set 
of dependent variables we choose becasue the result is unique up to the transformation (I3.18p . 

4. Averaged Operators 

In this section we collect some properties of the averaged operators Q and X>. 

4.1. Spectral Formulas for the Averaged Operators. Because A = Aq-Vx is skew-adjoint 
in H, its spectrum, Sp(^), is purely imaginary. Define 

cx{A) = {u eR : iu e Sp{A)} . 

For every u G o-{A) let denote the eigenspace associated with the eigenvalue iu of A and 
let be the H-orthogonal projection onto H^^. Then for every u G <y{A) one has 

{E^W, I W2)m = {W, I E^W2)m for every VTi, e H , 

E^ = E^ , = E^M . 

For every cu G a{A) let denote the smooth functions in — i.e. let = M^r] C°°(T'^). 
We will assume that each is dense in and is contained in the domains of V and Q. This 
assumption holds because we are working over the periodic domain T''. In general settings it 
would mean that every smooth eigenfunction of A must also satisfy any boundary conditions 
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associated with V, which is typically is not true. Our periodic setting we have explicit formulas 
for each E^^, however we will put off using them as long as possible. 
The Spectral Decomposition Theorem implies that 

e"*^ = e'^'^^E^ for every t G M . 

If G H has nonzero components only in a finite number of the M.^^ then there are a finite 
number of nonzero terms in the decomposition 

e-'^W = e~''^'E^W. 

If we assume moreover that each of these components is smooth then by using this decompo- 
sition in fl3.14p we obtain the spectral formulas 

Q{W, W)= Yl E^,+^,Q{E^,W, E^,W) , 
where we understand that E^-^^^^ = when Ui + 002 ^ (t{A). 

4.2. Bounds on the Averaged Operators. The spectral formulas (14. ip only apply to those 
Vl^ G H that have nonzero components in only a finite number of the EI;^, each of which is 
smooth. Denote this set by (SwC^. While this set is dense in H, it must be extended to larger 
classes of W. This is done by continuity once we obtain appropriate bounds on the forms 
associated with the operators Q and V. 

We begin by recalling some bounds on the forms associated with the operators A, Q, and 
V in terms of norms that are invariant under e"*'^. Specifically, we employ the H'^ spaces that 
are the completion of (Bi^C^ in the norms defined for every s G M by 



Here W{^) is the Fourier coefficient associated with the wave vector ^^7J^ . 

Lemma 4.1. There there exist positive constants Cjs, and Ct> such that for every Wi,W2 G 
C°°(T^) and every s' G M one has the bounds 

\{W,\AW2)J < C,,\\W^\\^^s'\\V,W2Us' . 

\{Wi\VW2)J < Ci,||V.W^i||e-.'||V.iy2||H^' . 

For every s > d/2 there exist a positive constant Cq such that for every Wi, W2, W3 G C°°(T'^) 
and every s' > one has the bounds 

\{W, I Q{W2,Ws))J < C'Q\\V,W^\\Ms\\W2\\m\\W3U, 

\{Wi I QiW2,Ws))J < r'CI^\\Wi\U.s'{\\V.W2\Us'\\Ws\\M^ + ||W^2||h=||V.W^3||h»') . 



(4.2) 



(4.3) 



Proof. We refer the reader to ^ for similar proofs. □ 
Because the norms are invariant under e"*-^ we can prove the following. 
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Lemma 4.2. For every s > d/2, s' > 0, and Wi, W2, W3 G ®ujC^ one has the bounds 

\{Wi\VW2)J < Cv\\V,W,\U..>\\V.W2\Us' , 
(4.4) \{Wi\Q{W2,W3))J <C^||V,M/^i||h^||W^2||m||W^3||h, 

\{Wi\Q{W2,Ws))J < 2^'C^||l^i||e-.'(||V.H^2||H='||W^3||H= + ||W^2||HH|V.H^3||e=' 
where Cv and Cq are the constants appearing in ( [^.^[ ) and of Lemma \4-1 ■ 



Proof. Because the M.^ norms are invariant under e and because e and commute, 
for every s > and every W G C°^{T'^) we have 

(4.5) ||e"*-^iy|| = \\W\\ , ||V,e"*-^PF|| = \\V,.W\\ . 

We first prove the bound on T> in (14.41) . From the definition of V given in (13.141) we see that 

{W, I VW2)^ = lim ^ r {W, I e'^Ve-'^W2)^dt 

J — 

= \im — [ (e-'^Wi\Ve-'-^W2)„dt. 
Hence, the bound on V in (14.21) and the invariances (14.51) imply 

< Jim ^ r \\V,e-'-^WA\JV,e-'-^W2\Ldt 



T^OO 2T J_rp 



r 

T^oo 2T 



hm ^ / ||V,W^i||H||V,H^2||edt = C^H V,W^i||e|| V,W^2| 

T 



This proves the bound on V in (14. 4p . 

We now prove the bounds on Q in (14. 4p . From the definition of Q given in (I3.14p we see that 

r-T 



hm ;1 r {e-'^W^\Q{e-'^W2,e-'^W^))^dt. 
— ^00 zi J 



im ^TT^ / 

T 

Hence, the first bound on Q in (14.31) and the invariances (14.51) imply 







1 


< 


lim 






T^oo 


2T 






ch 


< 


lim 








2T 










lim 








2T 



\Wo\um 



This proves the first bound on Q in (14. 4p . The second bound on Q in (14. 4p follows similarly 
from the second bound on Q in (14. 3 p and the invariances (14.51) . thereby proving Lemma W72[ □ 
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4.3. Averaged Quadratic Convection Operator. Until now the entropy has played no role 
in our analysis of the weakly nonlinear-dissipative approximation. It will now play a central 
role. Specifically, the fact system (12.11) has a strictly convex entropy implies that its flux F{U) 
has a special form [7l[T6]. We will use this form to show that the averaged quadratic convection 
operator satisfies a cyclic identity and some new bounds. 

4.3.1. Special Form of the Flux. Henceforth we will assume that the hyperbolic-parabolic sys- 
tem (12.11) has a strictly convex entropy that is thrice continuously differentiable over the convex 
domain U. Define the set V = {Hu{U) : U G U}. One can show that V C M'^ is a domain and 
that the mapping Hjj : Z// — > V is a diffeomorphism. For every G V we then define H*{V) G M 
and J*{V) G by 

(4.6) H*{V) = V^U - H{U) , J*{V) = V^F{U) - J{U) , 

where U eU is uniquely determined by Hu{U) = V. Because the entropy flux J{U) is related 
to H{U) and F{U) by ([53]), one can show that the mappings H* : V ^ R and J* : V ^ R'^ 
are continuously differentiable with 

(4.7) h;,{v) = u, ry{v) = F{u). 

It follows that F{U) is given by 

(4.8) F{U) = J^{Hu{U)) for every U eU. 

Because we have assumed that H is thrice continuously differentiable while F is twice contin- 
uously differentiable over U, it follows from (14.71) that both H* and J* are thrice continuously 
differentiable over V. This contrasts with J, which is twice continuously differentiable over U. 

4.3.2. Cyclic Identity. The fact that the flux F{U) of system (12. ip has the form (14.71) is central 
to our proof of the following identity. 

Lemma 4.3. (Cyclic Identity.) If system Ii2.1\) considered over the spatial domain T'^ has a 
strictly convex entropy H satisfying Iji2.4\ ) that is thrice continuously differentiable over U then 
for every Wi, W2, W3 G C°°(T'^) one has the cyclic identity 

(4.9) {Wi I Q{W2, W^3))e + {^2 1 Q{Ws, W,))^ + {Ws \ QiW,, W^))^ = , 
where Q is defined by Iji3.14\ )- 

Proof. We will first establish ( 14. 9 p for the particular Q associated with the conserved densities. 
In that case W = U. The extension of (14. 9 P to general Q then follows from the change of variable 
formula IKT8\f . 

We will begin by establishing (14.90 when Ui,U2,U^ G ©ujC'<2°- Once this is done, the extension 
of (14. 9 p to Ui, U2, f/3 G C°°(T'^) then follows by a density argument that uses the bounds on Q 
from Lemma [4. 2 [ 

Let Ui, U2, U3 G ®ojC^. From the definition of Q given by (I3.14p we see that 
(t/i I Q(f/2, f/3))e = ^1™^ ^ {Ui I e*-^Q(e-*-^f/2, e-'^Us))^dt 



(4.10) 



where f/i(t) = e~*-^[/i, U2^^ = e^*'^f/2, and U^it) = e^^-^U^ are quasiperiodic functions of t. 
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By (18121) and (13AD we see that 

{Ui\Q{U2,U3))^= [ UlGoV,-Qo{U2,Us)dx = - [ V.U^ ■ GoQo{U2,Us) dx , 

Jjd Jfd 

where Go = Huu{Uo) and Qo{U2, U3) = \Fuu{Uo){U2 ® U^). By differentiating (148|) twice and 
evaluating at Uo we obtain 

Fu{Uo) = Jvv{yo)Go. 

FuuiUo) {U2 ® ^^3) = Jvvviyo) {G0U2 ® GoU;) + Ji^^y (K>)i/c7f/C/(f/o) (^2 ® ^3) ■ 

where Vo = Hjj{Uo)- Hence, 

A = Ao-V^ = J^y(K)G'o-V,, 

and 

+ iV.^i^ ■ GoJ;.viyo)Huuu{Uo) {U2 ® f/s) 

= jvvviyo) {Go ■ ® G0U2 ® Got/s) + \HuuuiUo) [dA ®U2® U3) . 
Summing the above relation with its cyclic permutations gives 

V.U^ ■ GoQo{U2, U3) + V.Uj ■ GoQoiUs, Ui) + V.Uj ■ GoQo{Uu U2) 
= iV.- [j;.vviyo){GoUi ® G0U2 ® G0U3)] + \dt[Huuu{Uo){U^ ®U2® U^)] . 
Integrating this over T'^ we obtain 

{Ui I Q(^2, ^3))h + [U2 1 Q{U,, U,))^ + {U, I Q{U,, U2))^ 

at jjd 

The time average of this equation yields 

(f/l I Q(f/2, f/3))e + {U2 I Q(f/3, f/l))H + {Us I Q(f/l, ^2))h = , 

because the time average of the time derivative of a bounded function vanishes, thereby proving 
the lemma. □ 

The cyclic identity (14. 9 p yields the following bound. 

Lemma 4.4. If system \2.1\) considered over the spatial domain T'^ has a strictly convex entropy 
H satisfying ( [^.^[ ) that is thrice continuously differentiable overU then for every s > d/2 and 
every Wi,W2 e C°°(T^) one has the hound 

(4.11) \{Wi\Q{Wi,W2))^\ < \C'q\\W41\\S/.W2\\w . 

where Q{Wi, W2) is defined by \3.14^ and Gq is the constant appearing in lji4-3\ ) of Lemma \4.1 . 

Proof. The cyclic identity (SSI) and the fact that Q{Wi, W2) = QiW2, Wi) imply 

2{Wi I Q{Wu W2))^ + {W2 1 Q{Wi, Wi))^ = . 
It thereby follows from the second bound in (14.41) that 

\{W,\Q{W,,W2))J < '2\{W2\Q{W,,W,))J < '^G'Q\\W^U\\V,W2\\M^ , 
where Gq is the constant appearing in (14.31) of Lemma [4.11 □ 



12 N. JIANG AND C. D. LEVERMORE 

4.4. Averaged Dissipation Operator. Because the averaged dissipation operator T> is non- 
positive semidefinite, its Friedriclis extension is a self adjoint nonpositive semidefinite operator 
with domain Dom(X>). For every G H and every uo G (j{A) define to be the component of 
V in — so that 

V^ = E^Ven^, V= Y^V^. 

We claim that the nonpositve definite Hermitian form associated with V extends to the domain 
Herm(P) given by 

(4.12) Herm(^) = |y g H : - ^ (K, | VV^)^ < ooj . 

Recall that ^ Dom(P) for every uj G <y{A) and that T) is nonpositive over Dom(P), whereby 
(K; I 'PK;)h < for every u G a{A). Then for every V G Herm(X>) one has by (14.11) that 

(4.13) {V\VV)^= J2 (yu^lW^h < 0. 

We now characterize when the Hermitian form associated with V is nondegenerate. 
Lemma 4.5. 

(4.14) (y I < for every nonconstant V G Herm(P) , 
if and only if the Kawashima condition ^3. 7]) holds. 



Remark. This kind of characterization was first proved by Kawashima 

Proof. First we show that (14.141) implies the Kawashima condition (13. 7p holds. Suppose not. 
Then there exists a nonconstant V G H such that V G EIj^ for some uo G (y{A) and V G Null(P). 
But then V = V^e Herm(^) with 

{V I VV)^ = {V^ I WJh = {V I W)e = , 

which contradicts (14.141) . Therefore (14.141) implies the Kawashima condition (13. 7p holds. 

Next we show that the Kawashima condition ( 13. 7p implies ( 14.14p . Suppose that ( 14.14^ is 
false. Then there exists a nonconstant V G Herm('D) such that 

Q={V\VV)^= J2 (K.|WJh. 

However this is equivalent to 

(K. I T>V^)m = for every lu G cr{A) , 

which is equivalent to 

K; e Nul^P) for every u E a{A) , 

But V implies that K; 7^ for some u G a{A). But then for this u we have is 
nonconstant, ^ ^ui, and G Null(P), which contradicts the Kawashima condition (13. 7p . 
Therefore the Kawashima condition (13.71) implies (I4.14p . and the proof of Lemma is complete. 

□ 
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5. Strictly Dissipative Approximations 

The global existence theory presented in the next section will require that T) be strictly 
elliptic. When this is the case we say that the WND approximation is strictly dissipative. 
For many physical systems the WND approximation has this property. Indeed, the averaged 
dissipation operator has played an important role in the study of large time behavior of solu- 
tions to hyperbolic-parabolic systems of conservations laws and discrete Boltzmann equations. 
Kawashima [23, EE] showed that in one space dimension the averaged dissipation operator was 
strictly dissipative whenever the Kawashima condition fl3.7l) holds. In that case he showed that 
solutions of the original system are well approximated by solutions to an "effective artificial 
viscosity" system constructed using the averaged dissipation operator. Motivated by this idea, 
Hoff and Zumbrun [HI [19] studied multi-dimensional diffusion waves for the barotropic Navier- 
Stokes system through an artificial viscosity system constructed with the averaged dissipation 
operator. Recently, Bianchini-Hanouzet-Natalini [4] used the same idea to study the large-time 
behavior of smooth solutions for partially dissipative hyperbolic systems with a convex entropy. 
In each case the averaged dissipation operator was shown to be strictly dissipative through a 
detailed spectral analysis. 

To our knowledge there is no proof that the Kawashima condition fl3.7p implies that the 
averaged dissipation operator is strictly dissipative in multidimensional settings. Here we give 
a stronger criterion that does the job in our spatially periodic setting. In the Fourier represen- 
tation we have 

where ^o(0 ^"^^ ^oiO families of Go-symmetric matrices in M^^^ defined for every 

^ G M'^ by 

One then has 

whereby 

Our Kawashima-type criterion for strict dissipativity is given by the following. 
Theorem 5.1. Suppose that for some a > there exists f3 > such that 

(5.1) GoB£) + ^ AoiO^GoBoiOMO >/3Go for every f e . 
Then there exists 6 > such that 

(5.2) -{V\VV)^>6\\V.V\g for every V eC^iJ''). 

Remark. The Kawashima condition is satisfied whenever (15. ip holds. 
Our proof of Theorem 15.11 requires the following lemma. 

Lemma 5.1. Let G e C^^^ be Hermitian positive definite (i.e. G = G* > 0). Let A e C^^^ 
be G-Hermitian (i.e. GA = A*G) and B G C'^^'^'^^ be G-Hermitian nonnegative definite (i.e. 
GB = B*G > 0). Let Ga and Gb be constants such that \\A\\g < Ga and \\B\\g < Gb- Suppose 
that for some a > there exists (3 > such that 

(5.3) GB + i-A*GBA>pG. 
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Then there exists S > depending only on a, (3, Ca, o,nd Cb, such that 
(5.4) 1™;^/" e''^'GBe-'''^&t>5G. 

Proof. Let v G C'"^ be nonzero and define 

This is a smooth, nonnegative, quasiperiodic, real-valued function that will vanish whenever 
g-jtA^ is a null vector of -B. In order to prove (15. 4p we must obtain a lower bound for its 
average. The idea of the proof is to show that f{t) cannot be too small for long. 
The first two derivatives of f{t) are 

f{t) = iv*e''^\A*GB - GBA)e-''\, 

^^'^^ fit) = v*e''^* {2A*GBA - A*^GB - GBA^)e-''\ . 

Roughly speaking, we will show that f{t) is dominated by its first term when f[t) is small. For 
every rj > one has 

1 

The first term on the right-hand side above is just f{t)/rf while the second can be bounded by 
rfCXCBV*Gv. If we set rj^ = a^(3/{2CXCB) then we obtain the bound 

(5.6) \v*e''^\A*^GB + GBA^)e-''\\ < '^^^ f (t) + ^ v* Gv . 

a^p 2 

Let fi^ = {t G M : f{t) < ev*Gv}, where e > satisfies 
For every t G we obtain from (15. 3p . (15. 5p . (15. 6p . and (15. 7p the lower bound 



\v*e''^\A*^GB + GBA^)e-''\\ < —v*e'''^'GBe-''\ + 7]^v*e''^'A*^GBA^e'''\, 



f(t) > t;*e**^ {2a'GB + 2A*GBA)e-''\ - { 2a' + f{t) - -f v*Gv 

> a'(5v*Gv. 

For every t G i^'' we obtain from (15. 5p . (15. 6p . and (15. 7p the upper bound 

fit) < v*e''^' {2A*GBA)e~''\ + =^^4^ fit) + ^ v*Gv 

a'^p 2 



< {2ClCB + a'(3)v*Gv. 

Because / is continuous the set ^2" is open and is therefore a countable union of disjoint open 
intervals: 

^'=[j(albl). 

ken 
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Let (a, b) be any one of these intervals. Because f{t) satisfies the lower bound f lS.Sp while 
f(t) < ev*Gv for every t G (a, b), it is clear that the interval (a, b) must be bounded. 

We will begin by bounding b — a above and below. Because / is continuous over the bounded 
interval [a, b], and because a and b are not in fi", it follows that /(a) = f{b) = ev*Gv, and that 
/ takes its minimum at a point to G (a, b), at which /(to) = 0. Then 

(5.10) f{t) = f{Q+f{h)dh = f{Q+ / /(t2)dt2dti. 

to to 

We claim that 6 — a satisfies the bounds 

The upper bound of (15.111) is obtained from (I5.10p by using the lower bound (15.81) for f{t). 
For every t G (a, 6) we have 

fit)>fito) + l{t-tofa^Pv*Gv. 
Evaluating this at t = a and t = b yields 

ev*Gv > f{to) + |(to - afa^pv*Gv , 
ev*Gv > f{to) + 1(6 - toya^(3v*Gv . 

Because f{to) > while max{(to — a) ^, (6 — to) ^} > jib — a)'^, we obtain 

ev*Gv > \{b-afa^Pv*Gv, 

which yields the upper bound of (15.111) . 

The lower bound of (15.111) is obtained from (15.101) by using the upper bound (15.91) for f{t). 
For every t G (a, 6) we have 

f{t) < /(to) + |(t - Q\2GlGj, + 0^(3) v*Gv , 

Evaluating this at t = a and t = b yields 

ev*Gv < f{to) + |(to - ay{2ClCs + a^p) v*Gv , 
ev*Gv < /(to) + 1(6 - to)\2ClCs + v*Gv . 

Because min{(to — a)^, (6 — to)^} < j{b — a)^, we obtain 

ev*Gv < fito) + 1(6 - af{2GlCj^ + a^p) v*Gv , 

which yields the lower bound of (15. lip . 

Next, we bound the average of /(t) over (a, 6) from below. By again using the lower bound 
(ESD for /(t) in (15A0|1 we obtain 

' I" m it > m) + - + - - °> + - a^ff.-Gv 



b-aj^ 6 

> f{to) + l{b-afa^f3v*Gv. 

The lower bound of (15. lip then implies 
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We are now ready to prove (15.41) with 6 given by 



For every T > let = {k e N : (a^, bl) C (-T, T)} and 

k&K^ 

Then using the fact that f{t) > e over [— T, T] — fl"^, the lower bound (15.121) . and the fact that 
S given by (I5.13P satisfies 5 < e, we find 



2T 



7^1 f{t)dt> — 
1 

2T 



/ f{t)dt+ J2 ['fit)dt 



> 



iJ[-T,T]-n' 
meas([-T,T]-f]^)e+ (K ' ^1) 6 



v*Gv 



> 



meas ([-T,T] nfi^ -fi^) 



2T 



5v*Gv. 



The set [— T, T] fl fi^ — Vtfp is contained in the union of the (at most two) disjoint intervals 
(a^, h\) that contain — T and T. So its measure is bounded above by twice the upper bound for 
h — a given by (15.111) — namely, by 



meas([-T, T]r]Vt' - fi^) < 4 



2e 



Hence, we obtain 



2T 



7^ m dt > 



T 



2 

1 - - 



2e \ 2 



6v*Gv. 



Letting T — > cxd in this inequality yields (15. 4p . The limit of the left-hand side exists because 
fit) is quasiperiodic. □ 

We are now ready to prove Theorem 15.11 with the aid of Lemma 15.11 

Proof of Theorem 15.11 First observe that for every nonzero ^ G M*^ we set ^ = ^/l^l and 
have 

1 



-V{i) = lim 



T^OO 2T J_rp 

lim 



T^OO 2T J_rp 

|^|2 — /"^ e**^°(«)S„(Oe-'*^°(«) dt . 

T^oo 2T It 



Then by the fact Ao{C,) is Go-symmetric we have 



(5.14) 



-GoViO = leP lim — / G,e^*^°(«)5o(ae-^*^°«)dt 

T— >oo Z-Z 



lim 



1 



T^oo 2T 



WEAKLY NONLINEAR-DISSIPATIVE APPROXIMATIONS 17 

We now apply Lemma [5TT] with G = Go, A = Ao{^), B = Bo{^), 

GA = max{\\A£)\\G^ : f G S^-^} , = max { ||i?„(0||G„ : S^"'} • 

We find that for every ^ G S"'^^ we have the lower bound 

hm re^*^°(«)"G„5o(0e-^*^°(«)dt>5G., 

where S is given by (15 ■131) . Because 6 only depends on a, /3, Ga, and C^, it is independent of 
^. Combining this lower bound with (I5.14p yields 

-Go%{0>\^\^SGo for every e 
Then for every V G C^(T'^) the Plancherel identity implies 

-{v\vv)^ = - J2 vio*Gom)m>5 Yl \^\"vio*Gom = s\\v.v\\i. 

But this is (I5.2p . thereby proving the theorem. □ 



G 
lies 



6. Global Weak Solutions 
The weakly nonlinear-dissipative approximation is 

(6.1) dtW + AW +Q{W,W) = VW . 

Following the Leray theory for the incompressible Navier-Stokes system, we will show that if 
V is strictly dissipative then (16. ip has global weak solutions for all initial data W^^ G H. This 
result includes the Leray theory, so it cannot be improved easily. 

The key to obtaining global solutions in the Leray theory for the incompressible Navier- 
Stokes system is a so-called energy estimate. This designation is a bit misleading because, as 
we shall see, the estimate is better understood as an entropy estimate. 

6.1. Notion of Weak Solution. We call W G C([0, oo); w-H) nL^„^(dt; V) a Leray-type weak 
solution of weakly nonlinear-dissipative approximation (16. ip . if W satisfies the following weak 
form of (16. ip : 

= {V I W{h))^ - {V I W{h))^ - r {AV I W{t))^ dt 

(6.2) -^'^ 

- r ( ■GoQ{W{t),W{t))dxdt+ r f V,,V-GoB-V.W{t)dxdt, 

Jti JTd Jti Jfi 

for every function V G V. These solutions satisfy the entropy inequality 

(6.3) \\\W{t)\\l- !\w\VW)udt<\\\W'^\\l. 

Jo 

Of course, for every sufficiently nice W one has the identities 

(6.4) {W\AW)^ = 0, {W\Q{W,W))^ = 0. 
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6.2. Existence Theorem. The main theorem of this paper, the global existence of Leray-type 
weak solution to the weakly nonlinear-dissipative approximation (16.11) is as follows. 

Theorem 6.1. Let (16.11) be the weakly nonlinear-dissipative approximation about a constant 
state Uo of the hyperbolic-parabolic system h2.1\) with a strictly convex entropy ^2.4^ . Suppose 
that the diffusive operator D in (16. ip is strictly dissipative, i.e. that there exists 5 > such that 

(6.5) -{V\VV)^>6\\V.V\\l for every V eC\T''). 

Then, for every G V, there exists a solution W G C([0, oo); i(;-EI) fl Lf^^{dt;N) satisfying 
(IQ]) and fl63|l . 

Remark. If condition (15. ip is satisfied then V will satisfy (16. 5p and the above theorem insures 
the existence of at least one global weak solution. 

Remark. In Section U.2\ we apply our theory to the Navier-Stokes system of gas dynamics. 
The resulting averaged system includes the incompressible Navier-Stokes system as a subsystem. 
The question of uniqueness for system (16. ip is thereby at least as hard as that of uniqueness 
for weak solutions of the incompressible Navier-Stokes system. 

Proof. The strategy for our proof was introduced by Leray in the context of the incompressible 
Navier-Stokes system [28], see also [271 [36]. It is now classical compactness argument that has 
since been used to prove existence of global weak solutions for other equations [10] . It proceeds 
in four steps. We begin by constructing a sequence of approximate solutions. We then show 
that this sequences is relatively compact, first in some weak topologies and then in a strong 
topology. Finally, we show that limit points of this sequence satisfy (16. 2p and are thereby weak 
solutions of (16.11) . This strategy strikes a balance between the fact that compactness is easier to 
establish for weaker topologies and the fact that passing to the limit in nonlinear terms requires 
convergence in a strong topology. 

6.2.1. Step 1: Constructing Approximate Solutions. One can construct a sequence of approxi- 
mation solutions Wn by any method that yields a consistent weak formulation and an energy 
relation. Here we do this with the Galerkin method. 

Let {IHIn}^;^ be a sequence of subspaces of EI such that each H[„ lies within C°°{T'^), has 
dimension n, and satisfies EI„ C ]H[„+i. Assume moreover that this sequence is complete. Let 
Vn denote the orthogonal projection from EI onto IH[„. Completeness implies that for every 
y G EI one has VnV —>■ V a.s n —>■ oo. The Galerkin approximation of dimension n is the 
system 

(6.6) dtWn + VnAW^ + VnQiWn, Wn) = VaVW^ , 

where Wn takes values in IHI,„. This is a system of n ODEs. Its nonlinearities are quadratic, 
hence locally Lipschitz. The Picard existence theorem insures that system (16. 6p has local 
solutions. Taking inner product with PV„, and applying the identities in (16. 4p . we obtain the 
energy identity: 

(6.7) \\\Wn{t)\\l- !\wn\VWn)mAt' = \\\Wm, 

Jo 

for every t > 0. This energy identity immediately implies a global bound on the approximate 
solutions Wn, which thereby exists for all time. 
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6.2.2. Step 2: Establish Weak Compactness. We claim that the approximate solutions Wn are 
relatively compact in C([0, oo); w-H) fl w-Lf^^{dt; w-V). 

First, from the energy identity (16. 7p . we that Wn{t) is uniformly bounded in EI thus relatively 
compact in w-H for every t > 0. Next, from the dissipation property of V, see (16. 5p . we have 



(6.8) 



■ f WV^WJldt' < ~ f {Wn\VWn)^dt' <C. 

Jo Jo 



Thus Wn is relatively compact in w-Lf^J^ [0, oo) ; w-V) . We need only to verify that Wn is equicon- 
tinuous in C([0, cxd); w-H) thus by Arzela-Ascoli theorem the weak compactness is established. 
The equicontinuity can be derived from the weak form of the Galerkin system (16. 6p . 

[V I Wnit,) - Wnit,))^ = r {V I AWn{t))^dt 

Jti 

' {V I QiWnit), Wnit)))^dt + [\V\ VWnit))^dt , 



(6.9) 



for every function G V. We first prove the equicontinuity for test function V eYH C^(T'^), 
which is followed from the first and third bounds in the Lemma 14. 2[ Then we extend the class 
of test functions to V by standard density argument, thereby finishing the proof of Step 2. 

6.2.3. Step 3: Establish Strong Compactness. We claim that Wn is relatively compact in strong 
topology of Lf^^{[0, oo); H). It is a direct consequence of the weak compactness result in Step 2 
and the fact that the injection 

(6.10) C([0, oo); w-H) n w-Ll^{dt; w-V) — > Ll^{[0, oo); H) 

is continuous. 

6.2.4. Step 4-' Pass to the Limit. Step 2 ensures that there is a subsequence of Wn, which we also 
refer to as Wn, converges in C([0, oo); w-H) n w-Lf^^(dt; w-V) to a limit W E C{[0, oo); w-H) fl 
Lf^^{dt;Y). Step 3 ensures the convergence of Wn to W in Lf^^{[0,oo);'E.). All that remains is 
to show that the limit W satisfies the weak form (16. 2p as well as the energy inequahty (16. 3p . 
Toward this end we check convergence of each term in the respective regularized versions, (16.90 
and (16. 7p . respectively. Again, we first consider the test function V in the class Vfl C^(T"') 
then use density argument later. First 



[V I Wnit,) - Wnit,))^ ^ {V I Wit,) - W{t,))^ 



as n —y oo 



because of the relative compactness of Wn in C([0, oo); w-H). The convergence of the first term 
on the righthand side of (16. 9p is trivial. Note that 

\v\QiWn{t),Wn{t)))^dt=- r f V,V-GoQ{Wn{t),Wn{t))dxdt, 

Jti Jrd. 

and Q{Wn{t), Wn{t)) is quadratic in Wn- Thus the strong compactness of Wn in Lf^^{[0, oo); H) 
in Step 3 ensures the convergence of above term. We also note that 

[\v\VWn{t))^dt= r f V.V-GoB-V.Wn{t)dxdt, 
Jti Jti Jt<^ 

The convergence of above term is straightforward. Thus, we show that the limit W satisfy the 
weak form (16. 2p . thus is a weak solution to WND approximation (16.11) . 
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Now, to recover the energy inequality (16.31) from (16. 7p . first we note that for the initial data 
term 

The convergence of Wn in C([0, oo); w-H) and L^^^([0, oo); H), together with the fact that the 
norm of the weak limit is an eventual lower bound to the norms of the sequence, yields 

\\WU<\immf\\W^it)\g. 

n— >oo 

Similarly, the convergence of W„ in w-Lf^^{[0, oo); w-V) implies 

- f {W\ VW)m dt<- lim inf / {Wn \ VWn)m dt . 
Jo Jo 

Thus, we finish the proof of global Leray type weak solutions. □ 

6.3. Uniqueness Theorem. Uniqueness can never be asserted by such a compactness argu- 
ment, but generally requires the knowledge of additional regularity of the solution. For example, 
here we will prove the following weak-strong theorem. 

Theorem 6.2. Let Ui,U2 G C([0,oo); w-M) n Lf„,(dt; V) be two weak solutions of the WND 
system with initial dataU'^.U'^ G H. Let s > max{c^/2,l}. IfUi G L2([0, T]; )nLi([0, T]; V") 
for some T > then Ui G C([0,T]; w-M.^^^) and for every t G [0,T] one has the energy equality 

(6-11) wpml- f {uim,)^^i'=\\\u-\\i 

Jo 

and the stability bound 

(6.12) \\U,{t) - U^{t)\\^ < ^^v{c'q[ IIV.^i(t)IL.dt') \\U- - UT\\^. 

In particular, if = then U2{t) = Ui{t) for every t G [0,T]. 

Remark. Equation (16. lip is simply the assertion that the energy inequality satisfied by the 
strong solution Ui is in fact an equality. The bound (I6.12p is a basic weak-strong stability 
bound, from which the uniqueness assertion follows immediately. 

The key to the proof of Theorem 16.21 will be provided by the following lemma. 

Lemma 6.1. Let f/i, f/2 G C([0, oo); w-M.)r\Li^^{dt; V) be two weak solutions of the WND system 
with initial data U'^.U'^ G H. Let s > max{d/2,l}. If Ui G ^^([0, T]; ) n ^^([0, T]; V^) for 
some T > then Ui G C([0,T]; w-EI''"^) and for every t G [0, T] one has 



{U,{t)\U2{t))^ + / (Q(f/i,f/i)|f/2)jj+(t/l|Q(t/2,f/2))Hdt' 

(6.13) _ 

- / {^u,\u,)^ + {u,\vu2)^dt' = [uT\u'^)^. 

Jo 

This lemma will be proved later. Now we will use it to prove Theorem 16.21 

Proof of Theorem 16.21 The energy equality (16. lip follows by setting U2 = Ui in equation 
(I6.13P of Lemma Em using the cyclic identity (14. 9 p to see that (f/i | Q{Ui,Ui))^ = 0, and 
multiplying the result by |. 
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We now derive the stability bound f l6.12p . Add the energy inequahties for Ui and U2 and 
subtract equation fl6.13p from the result to obtain 



2 

(6.14) 



i||f/2(t) - f/iWlIn - r {QiUi, f/i) I t/2)e + {Ui I Q{U2, U2))^dt' 
Jo 

Jo 

Upon letting W = U2 - Ui (so that U2 = Ui + W), we see that 

(Q(f/i, U^) I U^)^ + (f/i I Q(f/2, f/2))e = (f/i I Q(f/i, f/i))e + {W \ Q(f/i, 
(6.15) + (f/i I Q(^i, ^i))jj + (f/i I Q(W^, f/i))jj 

+ {U,\Q{UuW))^ + {U,\Q{W,W))^. 
The cydic identity dS]) implies that (f/i | Q([/i, Ui))^ = and 

(W^ I QiUu U,))^ + (f/i I QiW, U,))^ + {U, I Q(t/i, l^))jj = . 
We thereby see that relation fl6.15l) reduces to 

(QiUu f/i) I f/2)e + {Ui I Q(f/2, f/2))e = (f^i I QiW, W))^ . 
When this relation is placed into fl6.14p . we obtain 

(6-16) ^llW^Wlle- f {U^\QiW,W))^dt' - [\w \VW)^dt' < IHW^^X, 

Jo Jo 

where W'"" = U'^ - U'"^. The third bound in (j43D gives 

|(f/i|QWW^))el<^al|V.f/i|L||W^||e- 
We combine this bound with the fact —{W \ T>W^^ > to see that (16.161) yields the inequality 

Jo 

The stability bound (16.121) then follows by the Gronwall Lemma. □ 

All that remains is to prove Lemma 16.11 

Proof of Lemma [mi We begin by showing that Ui G L2([0,T];H") n Li([0,r]; V") implies 
Ui G C([0, T]; w-H^^^). From equation (16.11) we see that for every W G H"'^'^"-'^) we have 

{W I U,{t2) - U^it,))^ = - r {W I AU,{t) + Q{Ui{t), U,{t)) - VU,{t))^dt . 

Jh 

We have the following estimates: 



t2 



{W\AU^it))^dt 



<c^1|i^IIh-(.-i) / \mt)\\usdt. 



t2 



tl 



By the first bound in (14. 3 p of Lemma 14.21 with s' = s we have 

-t2 



t2 
tl 



{W\VU,{t))^dt <Cv\\V,W\U^s \\V,Ui{t)\\usdt 

Jh 



<Cv\\W\U-,s^i, \\U,{t)hsdt. 
Jti 
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By the third bound in (14.31) of Lemma 14.21 with s' = s — 1 we have 



t2 



{W\QiU,it),U,it)))^dt 



<2^C-M|W^lle-(-) / \\Ui{t)\gsdt 



t2 



Because Ui G ^^([0, T]; H") n ^^([0, T]; V") it follows that {W \ Uiif))^ is a continuous function 
of t over [0,r] for every W G W'^ . Hence, Ui G C([0, T]; w-ff-^). 

We now prove that (I6J3D holds for every t G [0,T]. Let a G V{W^) and (3 G I'(]R) be 
moUifiers such that a > 0, /3 > 0, supp(/5) C (— cxd, 0] and 

/ a{x) dx = 1 , / pit) dt = l. 

For every e > define 6, eV{T'^x R) by 

For each i = 1,2 and e > define f/^e = ©e * t^i, so that for every x G T'' and t > we have 
Ui,{x,t) = {Q,*Ui){x,t) = / ee(x-x',t-t')f^i(a;',Odx'dt'. 
Each f/jf is a smooth function over M'^ x ]R_,_ that satisfies 

Because A is skew-adjoint we thereby see that 

= -{Qe*Q{UuU,)\U2e)^-{Uu\&e*Q{U2,U2))^ 

Upon integrating this equation over [0, t] we obtain 

{UUt) I f/2.(t))H + r (e, * t/i) I [/2.)e + {Uu I e, * Q{U2, t/2))e dt' 

^0 

- [\vUM.)^+{U,,\VU2.)^dt' 
Jo 

= {uuo)\U2M)u- 

We claim that 

lhn(f/i,(t)|f/2.(t))jj= (f/i(t)|f/2(t))H, 
lim(f/i,(0)|f/2.(0))j^= (f/r|[/-)jj, 

hm [\e.*Q{UuU^)\U2.)^dt'= f {Q{U,,U,)\U2)^dt\ 
''^^Jo Jo 

lim / {Uu\ee*Q{U2,U2))^dt'= [ {Ui\Q{U2,U2))^dt' , 
^^^Jo Jo 



(6.17) 



(6.18 
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Once these limits are established we can then pass to the limit in (16.171) to obtain (I6.13P and 
thereby complete the proof of Lemma [6TT1 The limits (16.181) are established by using the bounds 
from Lemma 14.21 and the convergence and boundedness properties of convolution. 

Now the first two limits are direct consequence of the convergence property of convolution in 
C{[0,T];w-M.^~^). Because Ui and U2 are continuous in time in w-EI*~^ and EI respectively, we 
have, for any t > 0, 

Uu{t)^Ui{t) in w-H^'-^) thus Uu{t)^Ui{t) in H, 
U2eit) ^ f/2(t) in w-H, 

as e — >^ 0, which imply that for every t > 0, 

(t/ie(t)|f/2.(t))e- (f/iW|f/2(t))e, 

as e 0. Thus we prove the first two limits. 
To prove the third limit, first, we have 



(6.19) 



(0e * Q{U^, U,) I f/s,)^ - (Q(t/i, t/i) I f/s)^ dt' 

/ {ee*Q{UuU,)-Q{UuU,)\U2e)^dt'+ [ {Q{UuU,)\U2,-U2)^dt'. 
Jo Jo 



Note that Uijs a strong solution, i.e. Ui E L°°([0,T];ff) n L'^{[0,T];Y') for s > d/2, and the 
structure of Q{Ui, Ui) is the derivative of Ui multiplying ?7i, thus 

Q(f/i,f/i) GL2([0,T];H). 

Then the convergence in i^^([0, T]; H) and the boundedness of mollifier imply that each integral 
on the right-hand side of (16.191) goes to as e — 0. 

We leave the proof of the fourth limit to the last step because it is the hardest one. We prove 
the fifth limit first. Applying the first inequality in Lemma 14.21 with s' = and the Holder 
inequality, we have 

I (VUu I f/2.)H - (^f^i I f^2)e dt' = /* (pUr, - VU, \ U2,)^ + {VU, \ U2, - f/s)^ dt' 
Jo Jo 

< Cv\\^x{Ule — t^l)||L2([0,T];H)||V2;f/2e||L2([o_r].e) 
+ C*©!! V2;t/i||L2([o,T] ;e)|| Vx(f/2e — f/2)||L2([0,r] ;e) • 

Note that V^Ui^ = (y^JJi)^, for i = 1, 2 , the boundedness of V^Uie in ( [0, T] ; H) , and the 
convergence 

{%Ui)^^VxUi in L2([0,T];H) as e , 

we finish the proof of the fifth limit. 

The main difficulty is that U2 is only a Leray weak solution, so Q{U2, U2) is not in L^, thus 
the method to prove the third limit is not applicable here. However, we have the following 
identity: for any functions U, V so that U^GoV e L\[0,T] ; L\dx)) and T > 0, 

rp rp 

(6.20) / [Q,*U\V)^dt= [ {u\e,*v)^dt. 

Jo Jo 
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Proof of (16.201) : By changing the order of integration, 

'(0,*f/|y)^dt 

/•T /• POO p 

= -^a{^)/3{'-f)U''ix',t')dx'dt'GoVix,t)dxdt 

Jo J J'' Jo Jf^ 

= r f r f ^a{^)p{'-^)V^{x,t)dxdtGoU{x',t')dx'dt' 

Jo J J'' Jt' Jjd 

= r {u\e,*v)^dt. 

Jo 

Applying fl6:2nD to Ui and Q(f/2, f/2), we have 

/ {Uu\Qe*Q{U2,U2))^dt' = [ {ee*Uu\Q{U2,U2))^dt'. 
Jo Jo 

Thus from the second inequahty in Lemma 14. 2^ 

/ {e,*Uu-Ui\Q{U2,U2))^dt' 

Jo 

C'qT^^ [ ||V,(e,*f/i,-f/i)|Uo^([0,T];ff)dt'||f/2||io 

Jo 



< 

Note that 



^(di;H) • 



v.(e,*f/,) = e,*(v.f/),. 

We claim that for any W G ^^([0, T]; H"), 

lim ||6e *We- W\\l2(io,t]-M') as e , 
which simply followed from the triangle inequality 

II ©e *We- W^||l2([o,t];H») < ||©e * " W^e || L2([0,T];H») + || W^e " ||L2([0,T];e'') , 

and the convergence property of convolution in L^([0, T]; H'^). We thereby prove the limit. 

We have now established all the limits asserted in fl6.18p and have thereby completed the 
proof of Lemma 16.11 □ 

7. Application to the Compressible Navier-Stokes System 

7.1. Compressible Navier-Stokes System. The compressible Navier-Stokes system of gas 
dynamics is an important example of a nonsingular hyperbolic-parabolic system with the strictly 
convex entropy. It governs the mass density p{x, t), bulk velocity u{x, t), and temperature 9{x, t) 
over i7 C M'' in the form 

dtp + V^- (pu) =0, 

(7.1) dtipu) + {pu®u + pl + S) =0, 
^t(|p|^r + P^) + "^x- {Ip\u\'^u + peu + pu + S ■u + q) = 0, 

where the specific energy e and pressure p are given by thermodynamic equations-of-state, 
e = e{p,9) and p = p{p,9), while the stress S and heat fiux q are given by the constitutive 
relations 

(7.2) S = -p{V:,u + {V^uf -^IV^-u) ~XIV^-u, q = -kVJ . 



WEAKLY NONLINEAR-DISSIPATIVE APPROXIMATIONS 25 

Here the coefficients of shear viscosity /x, bulk viscosity A and thermal conductivity k are given 
by formulas /i = ^{p,6), A = X{p,6), and k = k{p,6) that come either from a nonequilibrium 
(kinetic) theory or from fits to experimental data, while D is the dimension of the underlying 
microscopic world — usually D = 3. We require that D > max{2, d}. 

Equations (17. ip express the local conservation of mass, momentum, and energy. The con- 
stitutive relations (17. 2p for S and q model viscosity and thermal conductivity, which arise due 
to deviations of the gas from local thermodynamic equilibrium. Equations (17. ip reduce to the 
compressible Euler system when one sets S = and g = 0. 

The thermodynamic equations-of-state for the specific energy and pressure, e = e{p, 6) and 
p = p{p, 9), are assumed to be twice continuously differentiable over (p, 9) G and to satisfy 

(7.3) deeip, 6) > , dpp{p, 9) > , for every p > and ^ > . 
In addition, they must satisfy the Maxwell relation 

P%e + e% (0=0. 

This implies the existence of a function a = a{p, 9) that satisfies the differential relation 
This is equivalent to 

(7.4) da=id.--|^dp. 

We can identify a with the specific entropy. 

We make use of the convective form of system (17.11) . 

dtp + u- VxP + pVx ■ M = , 

(7.5) p{dtu + u ■ V^u) + V.^p + V^. ■ S = , 

p{dte + u ■ \7.£:) + pVa,. ■ M + S" : + Va,. ■ g = , 

and the differential specific entropy relation (17.41) to see that 

p p 11 

p{dtO + u ■ V^a) = -{dts + u- Vx-e) {dtp + u ■ V^.p) = --S -.V^u - -V^ - q . 

9 p9 9 9 

This can be put into the divergence form 

(7.6) dt{pa)+V^- (pMa + = -^S -.V^u - ^q-VJ . 

The local form of the second law of thermodynamics and the constitutive relations (17.20 imply 
that for any values of p, 6*, VxU + {Vxu)'^ — j^I^x - u, Vx - u, and Vx9 one has the inequality 

~S:VxU--^q- Vx9 = ^ I + {V^uf - :§/ V,- + A |V, - u\^ + k \Vx9\^ > . 

Because these values can be independently specified at any point in f2, the above inequality 
implies that p{p,9) > 0, X{p,9) > 0, and ii{p,9) > for every p > and 9 > 0. Of course. 
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these thermodynamic constraints are satisfied when jj, = X = k = 0, which is the case of the 
compressible Euler system. For the compressible Navier-Stokes system we require that 

(7.7) fi{p,e)>0, X{p,e)>0, K{p,e)>0, for every p > O and ^ > O, 

which are also consistent with the thermodynamic constraints. 

Now consider the system over a periodic box Q = T'^. Integrating the divergence form of the 
entropy equation (17.61) over T'^ yields 



(7.8) ^ / pa ax 



^ |V^M + (V^w)'^ - §J V^- + A |V^- + K \VJ\ 



dx . 



One sees from (17.71) and (17.81) that any stationary classical solution of the system for which 
p > and 9 > must satisfy 

Because u is periodic, one can use the first equation above to argue that V^u = 0. It then 
follows from (17.21) and (17.51) that Vxp{p, 9) = 0, which by the second equation above and (17.31) 
yields VxP = 0. This shows that the compressible Navier-Stokes system (l7.1H7.2p satisfies the 
nonsingularity condition (I2.10p . 

The compressible Navier-Stokes system (I7.m7.2l) is a nonsingular hyperbolic-parabolic system 
with a strictly convex entropy given by H{U) = —pa{p, 9) where U is related to p, u, and 9 by 

U = {p pui ■■■ puB IpIuI"^ + pe{p,9))^ . 

The set U is the range of this mapping restricted to the domain p > 0, u E W^, and 9 > 0. 
Whenever 9 ^ e{p, 9) is a strictly increasing function from onto M+ then U is given by 



I (f/o f/i ■ ■ ■ f/D Uu+if e M°+' : f/o > , 2f/of/D+i > U,^ + ■ ■ ■ + U^}. 



The function H : U M. will be strictly convex if and only if dge{p, 9) > and the sound 
speed is defined [17]. These conditions are satisfied by all thermodynamic equations-of-state 
that satisfy (17.31) . 

Finally, it is easily checked that the Navier-Stokes system (I7.m7.2p satisfies our Kawashima- 
type criterion (15. ip . whereby it also satisfies the Kawashima condition (13.70 . Moreover, because 
its weakly nonlinear-dissipative approximation will be strictly dissipative, our theory of global 
weak solutions applies to it. 

7.2. Weakly Compressible Navier-Stokes System. A fiuid dynamical system that for- 
mally includes both the acoustic and the Stokes systems is the so-called weakly compressible 
Stokes system 

dtp + Po^x ■ M = , 

Podtu + {{dpp)o p + {dgp)o 9) = po^x ■ [V^M + {^xuf - ^Vx ■ ul] + AoVx(Vx ■ u) , 
PoC^dt9 + 9o{dgp)o Vx-u = KoAJ , 

where = {806)0 is the specific heat capacity at constant volume. A so-called weakly com- 
pressible Navier-Stokes system that formally includes both the acoustic and the Navier-Stokes 
systems is the weakly nonlinear-dissipative approximation. It decomposes U into a component 
governed by the incompressible Navier-Stokes system and a component governed by a nonlocal 
quadratic acoustic equation that couples to the incompressible component. 
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(7.9) 



The incompressible component Um = {p,u,6) is governed by 

Po [dtu + u ■ V^u) + \7;p = Po^xU , 
PoC^{dt~e + u-Vj) =KoAj, 
where 

(7.10) V,-M = 0, id,p)op+idep)oe = 0, 

and denotes the specific heat capacity at constant pressure, while the acoustic component 
Uac = iv, V, x) is governed by 

(7.11) dtUac + AUac + V^-Qin{Urn, Uac) + V^'QaciUac, Uac) = VA^Uac , 

where 



(7.12) 



Vr A {; = , 



6o{dep)o ~ 



The diffusion coefficient v in fl7.12p is given by 
(7.13) 



2po 

where Cq is the sound speed defined by 

(7.14) cl = {d,p)o + 



2p^cY picyci 

Oo{dep)l 



Notice that z/ is positive if either po, Ao, or is positive! Q^.^ and Qac are nonlocal operators 
defined in the following way. We first introduce the orthonomal basis of the acoustic mode 
Null(^)^: 

^ \ 

(7.15) H^{x) = ^ 



^\k\ 



Ak ■ X 



It is easy to check that is the eigenvector of the acoustic operator A with eigenvalues 
±iCo|fc|. Then any acoustic component Uac can be represented as 

(7.16) Uac = J2UkHH^). 



where U^ is the coefficient of Uac with respect to the basis under the inner product 
The nonlocal operator Q^^ can be written as 



(7.17) 
where 
(7.18) 



QiniUin,Uac) — ^ A^(t/j„)iJj^ 
S,m 



A^(f/i, 



itfciti=m 
|A;| = |m| 



^^{urm)k + {k.m)u^ + ^(c^k + c,m) 



\m\ 



\m\ 
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The nonlocal operator Q^^ can be written as 

(7.19) Q.n(^ac,^ac)=E^-^-(^)' 

m 

where 

(7.20) xi = c, Y: u^u^iUrr 

k+l=m 
±{k)\k\+±{l)\l\=±\m\ 

where the constants ci, C2, C3, and C4 are calculated in p2]. 

We have shown that the weakly compressible Navier-Stokes system has global weak solutions 
in L^. This result includes the Leray theory, so it cannot be improved easily. As with the 
Leray theory, the key to this result is an "energy" dissipation estimate. We have, (extending 
ideas of Masmoudi and Danchin) also adapted a Littlewood-Payly decomposition to show the 
acoustic part is unique for a given incompressible component [21j. Moreover, we have derived 
the weakly compressible Navier-Stokes system directly from the Boltzmann equation for the 
case p = p9 and e = ^9 [22] . 
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